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We have performed elastic diffuse neutron scattering studies on the relaxor Pb(Mgi/3Nb2/3)03 
(PMN). The measured intensity distribution near a (100) Bragg peak in the (hkO) scattering plane 
assumes the shape of a butterfly with extended intensity in the (110) and (lIO) directions. The 
temperature dependence of the diffuse scattering shows that both the size of the polar nanoregions 
(PNR) and the integrated diffuse intensity increase with cooling even for temperatures below the 
Curie temperature Tc ~ 213 K. 

PACS numbers: 77.80.-e, 77.84.Dy, 61.12.Ex 



I. INTRODUCTION 



II. EXPERIMENT 



o 

> 
O 

00 

o 
m 
o 



I 

o 
o 



X 
S3 



Pb(Mgi/3Nb2/3)03 (PMN) is a typical relaxor fer- 
roelectric that has a broad and strongly frequency- 
dependent dielectric constant. This system has attracted 
much attention in the last few years, because of the dis- 
covery of an ultrahigh piezoelectric response in solid so- 
lutions with PbTiOa (PT) near the morphotropic phase 
boundary (MPB)i. Pure PMN is considered to be a pro- 
totype relaxor, and has been studied extensively in recent 
years^. Nevertheless, certain aspects of this system are 
not fully understood. Unlike one of its close analogues, 
Pb(Zni/3Nb2/3)03 (PZN), PMN remains cubic below the 
(field-induced) Curie temperature Tc ~ 213 K'^-'*, and 
was believed to exhibit no macroscopic ferroelectric phase 
transition in zero field. 



Efforts to understand the dynamical properties of re- 
laxors have focused on the polar nanoregions (PNR) 
present in these compounds. Burns and Dacol^ found 
that on cooling such systems, PNR start to form at 
temperatures a few hundred degrees above Tq- This 
temperature was later called the "Burns Temperature" 
Td, Td ~ 620 K for PMN. Neutron scattering measure- 
ments®'^ have demonstrated that diffuse scattering ap- 
pears in PMN between 600 K and 650 K, consistent with 
previous optical measurements by Burns and Dacol. The 
diffuse intensity that develops in relaxors below Td has 
been identified with the PNR. Various neutron and x-ray 
measurements on diffuse scattering have been carried out 
in order to investigate how PNR are formed, and to de- 
termine average sizes and polarizations (atomic shifts) 
at different temperatures^-^-^"-^^. Phonon contributions 
have always been a potential source of contamination 
in diffuse scattering measurements. In this paper, we 
present, and attempt to interpret, elastic diffuse scatter- 
ing measurements on pure PMN. 



The experiment was carried out using the Disk Chop- 
per Spectrometer (DCS)^^ at the NIST Center for Neu- 
tron Research. The sample was a high quality single crys- 
tal of PMN with ^ 0.1° mosaic, and a mass of 3.25 g, 
grown at the Simon Fraser University in Canada. A num- 
ber of neutron scattering studies using this crystal have 
already been published^ii^ii^ii^. The room temperature 
lattice parameter is a = 4.04 A. At each of several tem- 
peratures, time-of-flight spectra were collected for each of 
at least 51 successive crystal orientations 0.5° apart us- 
ing 325 detectors with active dimensions in and normal 
to the scattering plane of 31 mm and 400 mm respec- 
tively. The DCS detectors are located 4000 mm from the 
sample at scattering angles from 5° to 140°. At an inci- 
dent neutron wavelength of 5.5 A, the FWHM (full width 
at half maximum) of the elastic resolution function was 
^ 0.085meV, and the time between pulses at the sample 
is 9 ms, effectively eliminating frame overlap. Data were 
collected near a (100) Bragg peak in the (hkO) scattering 
plane. The choice of a low Q refiection helps minimize 
phonon and multiphonon contributions. 



III. RESULTS AND DISCUSSION 

In Fig. n we show the neutron diffuse scattering 
around a (100) peak in the (hkO) plane measured at dif- 
ferent temperatures. Compared with similar x-ray mea- 
surements, our neutron scattering measurements have 
much higher energy resolution, and almost all phonon 
contributions can be removed by integrating the scat- 
tering over an appropriate energy range, 0.075 meV ei- 
ther side of the elastic peak position. In fact, the energy 
spectra show that there is no significant scattering in- 
tensity outside this integral range, i.e., phonon contribu- 
tions are indeed very small around the (100) Bragg peak. 
The remaining low energy phonon contributions included 
within this energy integral (phonons with energy trans- 
fer \fiio\ < 0.075 meV) can also be estimated, because 
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FIG. 1: Logarithmic plots of the neutron elastic diffuse scat- 
tering intensity around a (100) Bragg peak, at different tem- 
peratures. 

phonon intensity increases with temperature, while the 
elastic diffuse scattering intensity decreases. Since no 
change in lineshape with temperature was observed, we 
believe that phonon contributions to the elastic diffuse 
scattering intensity are negligible. 

Comparing the data sets plotted in Fig. ^ one can 
clearly see that diffuse scattering develops around the 
(100) peak with decreasing temperature. At ^ 650 K, 
the "butterfly" -shaped diffuse scattering intensity pat- 
tern observed at low temperatures is absent, and only 
a small trace of intensity transverse to the wave vec- 
tor Q = (1,0,0) remains. An enlarged plot of the dif- 
fuse scattering intensity measured at 200 K is given in 
Fig. 121 The diffuse scattering intensity has a "butterfly" 
shape and extends from the Bragg peak along (110) and 
(110) directions. Previous x-ray diffuse scattering mea- 
surementsi^ show that the polarizations of the PNR are 
along {111} directions in PMN, and the local symmetry 
is rhombohedral. This {111} type polarization gives the 
"butterfly" shaped diffuse intensity pattern in the (hkO) 
plane. The present elastic diffuse scattering results are in 
good agreement with the x-ray results. Recent neutron 
diffuse scattering measurements in the (hhk) scattering 
plane by S.-H. Lee et alr '^ provide additional evidence to 
support the {111} type polarization in PMN. 

Fig.|3|shows cuts through the data along the (110) di- 
rection at three selected temperatures. At 400 K, which 
is well above Tc, the diffuse scattering intensity is al- 
ready measurable. At higher temperatures, the diffuse 
scattering intensity becomes even weaker and harder to 
measure. We use a simple Lorentzian function to describe 
the spatial correlation of the atomic displacements con- 
tributing to the diffuse scattering: 

" 4q2 + r2)' 

where Iq is the integrated diffuse scattering intensity, P = 
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FIG. 2: A smoothed logarithmic plot of the neutron elastic 
diffuse scattering intensity at 200 K. 

1/^ is the inverse of the real space correlation length ^, 
and q is the length of the wave vector measured from the 
(100) Bragg position. Fig. |31 shows the results of fits to 
the Lorentzian, plus a Gaussian that describes the central 
Bragg peak, and a flat background. We have been able to 
obtain good fits using this model. The fitting parameters 
are the integrated intensity Iq and the half width at half 
maximum P of the Lorentzian, the intensity and width 
of the central Gaussian, and a flat background. Iq and 
^ = 1/P are plotted in Fig. El 

In Table, m we list the integrated diffuse scattering in- 
tensity Jo, correlation length ^, and the intensity at a 
selected wave vector Q — (1.05,0.05,0) vs temperature 
T. The elastic diffuse scattering starts to be noticeable 
around T -400 K, much below Td « 620 K. The "corre- 
lation length" ^ is a direct measure of the length scale of 
the static PNR. According to our results, the PNR are 
small when they first appear at high temperatures, with 
average sizes around 15 A (see Fig. 0J. Both the spatial 
correlation length of the atomic displacements and the 
integrated intensity of the diffuse scattering increase on 
cooling, even at temperatures below Tc (Fig. At low 
temperatures the length scale of the PNR reaches '^65 A. 

From T = 300 K to 100 K, the volume of a single PNR 
{V cx increases by a factor of ~60, yet the integrated 
intensity only increases by a factor of ~10. Writing Jq as 
the product A^^'^IQ • where N is the total number of 
PNR and d is the average displacement of atoms within 
the PNR, we conclude that A^|Q-5p increases on cooling 
from high temperatures and then drops dramatically at 
around Tc, remaining roughly constant below Tc- This 
is illustrated in Fig. \E\ which shows Iq /S,^ as a function 
of T. A likely scenario is that on cooling below Tc, 5 
remains relatively constant whereas N decreases as the 
smaller PNR merge together. 

Previous inelastic neutron scattering measurements on 
PMN^"* have shown that the transverse optic (TO) mode 
phonon become overdamped near the zone center, start- 
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TABLE I: Integrated diffuse scattering intensity 7o, correlation length ^, and the intensity at a selected wave vector Q = 
(1.05, 0.05, 0) vs temperature T. 



T (K) 


10 


100 


150 


200 


250 


300 


400 


500 


600 


Integrated Intensity Jo 


504.9 


437.6 


361.4 


287.1 


145.0 


59.3 


5.62 


3.06 




C(A) 


62.4 


64.5 


53.8 


47.5 


19.3 


14.6 


11.2 


13.3 




-^Q=(l. 05, 0.05,0) 


174.6 


163.7 


149.6 


137.0 


130.2 


87.22 


20.63 


7.90 


2.03 
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FIG. 3: Diffuse scattering intensities along the (110) direc- 
tion near the (100) Bragg peak, measured at 10 K, 200 K, 
and 400 K. The intensity profiles are fitted using a broad 
Lorentzian function (solid lines), a narrow Gaussian (dashed 
lines) that describes the resolution-broadened Bragg peak, 
and a flat background. 



ing at Td ~ 600 K, which is roughly the temperature at 
which the PNR start to appear. We believe that the over- 
damping of the soft TO mode phonon is directly associ- 
ated with the formation of the PNR. In other words, the 
PNR may originate from the condensation of the soft TO 
mode. Reexamining the neutron diffuse data on PMN by 
Vakhrushev et ali^, Hirota et al& have proposed a new 
interpretation of the atomic displacements derived from 
the data. The atomic displacements can be decomposed 
into the sum of two terms with comparable magnitudes: 
Ssoft -'r 5 shift- Ssoft values satisfy the center-of-mass con- 
dition and are consistent with the values derived from 
inelastic scattering intensities from the soft TO mode. 
The other term, S shift, shows that the PNR are shifted 
along their polar directions relative to the surrounding 
lattices. This "uniform phase shift" has established con- 
vincing connections between the diffuse scattering from 
the PNR and the condensation of soft TO mode phonons. 
It is important to note that the magnitude of 5 shift is 



FIG. 4: Top frame: the correlation length ^ as a function of 
T. Bottom frame: integrated intensity Iq as a function of T. 



comparable to Ssoft , which can be estimated to be around 
1/10 of the lattice spacing. This large phase shift creates 
a huge energy barrier, preventing these polar regions from 
merging into the surrounding lattices. 

Recently, Gehring et al. have performed measurements 
that unambiguously established a true ferroelectric soft 
mode in PMN^'^. Further inelastic neutron scattering 
studies by Wakimoto et ali^ show that the soft mode 
extends to much lower temperatures below Tc, and that 
the phonon energy squared (ftwg)^, which is inversely pro- 
portional to the dielectric constant e, increases linearly 
with decreasing temperature. These results represent a 
dynamical signature of a ferroelectric phase below Tc- 
The phase below Tc is now believed to be the recently 
discovered phase XiSiSS, which exhibits an average cu- 
bic structure. However, unlike conventional ferroelectric 
phase transitions, no macroscopic (rhombohedral) lattice 
distortion has been observed. The decoupling of the lat- 
tice distortion from the ferroelectric polarization is quite 
intriguing, and we believe that the interaction between 
the PNR and the surrounding lattices is the key in solv- 
ing this puzzle. 

Our current understanding of PNR in relaxor systems 
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FIG. 5: Plot of 7o/?^ vs T. Parameters were obtained from 
the fits to the diffuse data. 



can be described as follows: At the Burns temperature 
Td, the soft TO phonon mode becomes overdamped near 
the zone center, and starts to condense into PNR. These 
regions form with local polarizations along {111} direc- 
tions and are shifted uniformly along their individual po- 
larization direction. The number and size of PNR in the 
system increases with cooling. At the phase transition 
T = Tc, a large scale overall "freezing" of the PNR oc- 
curs. Small PNR merge into larger ones and the total 
volume of PNR in the system keeps increasing. The re- 
lated ferroelectric soft mode lifetime increases below Tc, 



and the overdamping near the zone center disappears. 
A macroscopic ferroelectric polar phase without lattice 
distortion is then established. Below Tc the size of the 
PNR can grow slowly with further cooling. However, if 
the coupling between the PNR and the surrounding lat- 
tice is not sufficiently strong, as is the case in pure PMN, 
then the energy barrier created by the uniform phase shift 
would prevent the PNR from merging further and form- 
ing macroscopic lattice distortions. The resulting phase 
will have a polar lattice of average cubic structure, but 
with embedded (rhombohedrally) polarized PNR. 

In addition to neutron and x-ray diffuse scattering 
measurements, recent Raman studiesSi and specific heat 
measurements^ have also provided useful information on 
PNR in PMN. We are planning to carry out further dif- 
fuse scattering measurements on a series of relaxor sys- 
tems in the near future. 

In summary, we have shown high resolution neutron 
elastic diffuse scattering data from PMN in the tempera- 
ture range of 10 K to 650 K. The formation and develop- 
ment of PNR at temperatures above, around, and below 
Tc have been carefully studied. A merging of smaller 
PNR into larger ones occurs at the ferroelectric phase 
transition, while further growth of the PNR into macro- 
scopic rhombohedral domains can not be achieved below 
Tc. 
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